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Summary 

A method is described for designing an aero- 
dynamic configuration for a specified performance 
vector, based on results from several similar but not 
identical trial configurations, each defined by a geom- 
etry parameter vector. The theory indicates that the 
method is effective provided that (1) the results for 
the trial configurations provide sufficient variation so 
that a linear combination of them approximates the 
specified performance and (2) the differences between 
the performance vectors (including the specified per- 
formance) are sufficiently small that the linearity as- 
sumption of sensitivity analysis applies to the differ- 
ences. A computed example describes the design of 
a high supersonic Mach number missile wing-body 
configuration based on results from a set of four trial 
configurations. 

Introduction 

Following the discussion in references 1 and 2 of 
the problems inherent in most design-by-optimization 
procedures, the method described in reference 3 was 
developed to avoid most of these problems. The 
method of reference 3 was developed specifically for 
application to problems involving large (in running 
time and/or storage) analysis codes so that a very 
limited number of analysis runs are practical. It at- 
tempts to take maximum advantage of the designer’s 
experience and intuition. Although the problems dis- 
cussed relate to aerodynamic design, the method is 
applicable in a more general context. 

The method described in reference 3 is specifically 
applicable to problems in which a shape function is 
designed to obtain a specified pressure distribution 
corresponding to that shape. For many aerodynamic 
problems, especially for synthesized configurations, 
the shape functions are not specified directly but by 
a set of shape geometry parameters. The desired per- 
formance may not be specified by the pressure distri- 
bution as such, but by integrated quantities such as 
lift, thrust, moments. Furthermore, it is important 
not to generate a sensitive “operating-point” design 
or to generate a design based on optimizing a sin- 
gle payoff quantity, such as lift-drag ratio for a single 
flight condition without consideration of off-design 
conditions. It is preferable to specify a reasonable 
goal in terms of performance parameters designated 
over some range of flight conditions (i.e., a perfor- 
mance index). This type of design problem is treated 
in the present analysis. 

The method to be described does not generate an 
optimum configuration in an absolute mathematical 
sense. However, when the basic assumptions of the 


theory are met, it makes optimal (in a true mathe- 
matical sense) use of the available computed results 
to generate an improved design. It does not require 
the expensive calculation of sensitivities for the vari- 
ation of each individual geometric shape parameter. 
Finally, it avoids the pitfalls of single-point-design 
configurations by specifying the goal to be obtained 
in terms of a performance index, or vector, rather 
than by a single performance parameter. 

Symbols 

a^di coefficients 

b{ coefficients determined by equa- 

tion (10) 

Cj coefficients or parameters 

Cd drag coefficient based on inviscid 

drag only 

Ci lift coefficient 

d vector denoting difference between 

two performance parameter vectors 

set of orthonormal vectors 

g vector whose components are an or- 

dered set of parameters prescribing 
configuration geometry 

Hi , hi coefficients of polynomial terms 

L/D = C L /C D 

li lengths specified in inlet example 

(fig- 1) 

m number of vectors 

M free-stream Mach number 

n number of components in each 

vector 

p performance vector 

P predicted (in figures) 

S specified (in figures) 

t normalized geometry variable 

w weight factor 

z, j/, 0 Cartesian coordinates 

a angle of attack 

A, A functions of z determining size 

variation 

a function of x determining shape 

variation 


Subscripts: 

d desired or design 

i, j. k indices 

Analysis 

The term configuration is used in this report to 
refer not only to a complete vehicle but also to any 
component whose performance is calculated from its 
geometry and the existing flight or flow conditions. 
It is assumed that small changes in the configuration 
geometry result in small changes in performance. 
This is the basic assumption of sensitivity analysis. 
It is noteworthy that this assumption holds even in 
the presence of shock waves if the shock locations 
and strengths change slightly with small geometry 
changes. 

The method can be described as a step-by-step 
procedure. The basic idea of each of the six steps is 
as follows: 

1. Compute the performance quantities for sev- 
eral configurations, each defined by assigning values 
to a set of geometry parameters. Compute weights 
for the performance quantities to render them com- 
mensurable (comparable in magnitude). The result- 
ing performance vectors for these base configurations 
are termed the “basic vectors.” 

2. Establish a desired performance vector as a 
goal. 

3. Approximate the desired performance vector 
formally as a linear combination of the basic perfor- 
mance vectors. Express this relationship in terms of 
incremental vectors. 

4. Determine the coefficients in this relationship 
so that the approximation is optimal. 

5. Use these coefficients to determine a set of 
geometry parameters representing the optimal com- 
bination of the original configurations in accor- 
dance with the local linearity assumption of sen- 
sitivity analysis. Then analyze this synthesized 
configuration. 

6. Compare the results of this analysis with the 
desired results. Then, if necessary, add more config- 
urations to the base, and repeat the procedure. 

These steps are now discussed in detail. 

Step 1: The first step in the procedure is to ana- 
lyze several practical, and hopefully good, configura- 
tions which possess a degree of geometric similarity. 
For example, it is not feasible to compare a forward- 
swept-wing fighter with several swept-back configu- 
rations with large sweep angles. Each configuration 
must satisfy all the required geometric constraints. 

The geometry for each configuration is defined by 
a set of shape parameters which is ordered to form a 


geometry vector g. The corresponding performance 
parameters (lift, drag, moments, etc., at various 
operating conditions) are ordered into an array p, 
which is denoted the performance vector. 

There are two aspects to weighting performance 
parameters: importance and magnitude. The pos- 
sibility of weighting the performance quantities so 
as to deliberately bias the result in accordance with 
the relative importance of the various performance 
criteria is a highly problem-dependent consideration 
and is to a certain extent subjective, and therefore is 
not considered here. However, the second aspect of 
weighting must be considered in multiobjective op- 
timization problems since the various performance 
parameters are, in general, not commensurable (or 
comparable); that is, they may be of different or- 
ders of magnitude or of different dimension. Thus, 
the smaller quantities have less influence in the op- 
timization procedure, and the final design is biased 
toward the larger performance criteria, even though 
the smaller quantities are of equal importance. An 
obvious way to deal with this problem is to weight the 
various components of the performance vector. The 
weights can be determined as follows: (1) compute 
the root-mean-square (rms) value of the ith compo- 
nents for all configurations, (2) form an average vec- 
tor from these rms values, (3) normalize this vector, 
and (4) take the reciprocals of the zth component 
as the corresponding weight to be applied to the ith 
component of each performance vector. Thus, for 
example, if there are m vectors each with n compo- 
nents, the weight factor for the kth component may 
be computed by rms averaging. Denoting the ith 
component of the jth vector by p t y, 


w k 


n m 

E E v% 

t=ij=i 

m 

n £ plj 

j-l 


(i) 


Henceforth, it will be assumed that the performance 
vector p consists of weighted components. 

Step 2: The next step in the procedure is to es- 
tablish a goal in the form of a desired performance 
vector which represents an improvement over the 
best of the configurations analyzed. If the configura- 
tions which have already been analyzed were actually 
“good” configurations, then a modest improvement 
over the best of these should represent a satisfactory 
goal. In this case, one proceeds to the next step in 
the procedure. 

Suppose, however, that the actual required per- 
formance is well beyond that of any of the config- 
urations analyzed. Then, one may proceed by one 
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of several alternatives. If one can surmise from the 
results of the initial calculations the nature of the 
defect in the configurations analyzed, this knowledge 
may lead to the scrapping of the original type of con- 
figuration and starting anew with a different class of 
geometries. 

If, on the other hand, a study of the results of the 
initial calculation indicates a way in which the basic 
geometry type can be modified to gain large steps 
toward the performance goal, then these improved 
versions can be added to the base of configurations 
to be analyzed. 

Finally, if the calculations provide no obvious clue 
toward making a significant design improvement, one 
can work gradually toward the goal by first specifying 
a subgoal representing say 20 percent of the difference 
between the best of the configurations and the actual 
required goal. Then after taking this initial step, 
add more configurations to the base, specify the 
next subgoal, and continue in this manner until the 
required goal is achieved. 

Step 3: The next step in the procedure is to 
approximate the desired performance vector by a 
linear combination of the performance vectors of the 
basic configurations 

m 

5 iPi ( 2 ) 

1 

subject to the constraint 

X> = 1 (3) 

1 

where the coefficients are still to be determined. This 
constraint preserves the scaling; that is, it assures 
that the components of p d will be close in magnitude 
to the corresponding components of p t , as is required 
for linearity in differences. The significance of this 
constraint in terms of the geometry parameters is 
discussed subsequently in step 5. 

Equation (3) reduces the number of required co- 
efficients in equation (2) by one, so that equation (2) 
can be written in terms of n — 1 coefficients a i by 
setting 


71—1 

Q,n = 1 ^ ^ 

1 


(4a) 

o>i cl^ (i 1, 2, . . • , 

n — 1) 

(4b) 

Thus, equation (2) becomes 



Pd « a ^j ? n + 

n-1 

^ a iPi 
1 

(5) 


In this form, it is easily seen that equation (3) is 
satisfied. Equation (5) can also be written in the 
form 

n- 1 

Pd^Pn+Yl a i(P n ~ Pi) ( 6a ) 

1 

or 

ra-l 

Pd~Pn^J2 aiipn - Pi) (6b) 

1 

In this form, p n is sometimes called a basis vector 
while the remaining vectors p{ (i — 1, 2, . . . , n — 1) 
are called comparison vectors. • 

Equation (6b) is an equation for the approximate 
representation of the difference p d — p n and can be 
written, with difference notation, as 

71-1 

dd = £ Mi (7) 

1 

Now since the input geometry parameter vectors 
gi define configurations of the same basic type that 
each satisfies all the prescribed geometric constraints, 
then the differences g n — g i (i = 1, 2, . . . , n — 1) 
should be small (have small components). Further- 
more, since each of the configurations analyzed is 
assumed to represent a good design, whose perfor- 
mance vector roughly approximates p d , then the dif- 
ference vectors d d and d z (i — 1 , 2, . . . , n-1) should 
be small. These are simply statements of some of the 
sensitivity analysis requirements. 

Step 4: It yet remains to determine the coeffi- 
cients in equation (7). An optimal choice of these 
coefficients would provide the closest approximation 
(in the least-squares sense) to d d that is possible 
with the n-1 vectors d z . If the vectors d; were an 
orthonormal set, optimal coefficients would be ob- 
tained simply by projecting d d onto each d{. In gen- 
eral, the vectors d z are not orthonormal, but they 
can be orthogonalized by the Gram-Schmidt process 
and normalized. That is, a new set of n - 1 vectors 
e; can be defined as a linear combination of df 

i 

~ ^2 c ijdj (^) 

3= 1 

such that the e ^ vectors form an orthonormal set. 
The procedure for computing the coefficients c^- 
in equation (8) is well-known but is given in the 
appendix in the present notation for the convenience 
of the reader. Now d d can be optimally approximated 
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by a linear combination of the orthonormal vectors 

n— 1 

dd, ^ (9) 

l 

with the coefficients 6; determined by 

k = d d • e { (10) 

With the coefficients b i determined in this manner, 
the right-hand side of equation (9) represents the pro- 
jection of dfi onto the n - 1 dimensional linear sub- 
space determined by the vectors This projection 
defines the closest point of the surface to d Since 
the same subspace is determined by the vectors d*, 
this closest point can also be represented as a linear 
combination of the vectors d z . This is accomplished 
by substituting from equation (8) into equation (9): 

n-1 

d d ^ y 1 b{€ 2 
1 

n— 1 i 

— b{ C ikdjc 
1 1 

Tt — 1 n— 1 

= ^ ^ c ik^k i c ik — 9 for k > i) 

1 1 

n— 1 /n — 1 \ 

= ( E b i c ik ) d k (c lk = 0 for k> i ) 

k= 1 V z!= 1 / 

71 — 1 

= E a k d k (11) 

l 

with 

k 

dk = Y^ b i c ik ( c ik = 0 for k > i ) (12) 

7=1 

These are the optimal coefficients to be used in 
equation (6b). 

Step 5: Now, with these same coefficients we syn- 
thesize a configuration geometry from the geometry 
vectors gf 

71-1 

9d = 9n~J2 a i^ n ~ ( 13 ) 

1 

With the original assumption that the differences and 
the sum of the differences are small, it follows from 
the linearity assumption of sensitivity analysis that 


the performance of this configuration satisfies equa- 
tion (6a) within the accuracy of the local linearity 
assumption. 

Equation (13) can also be written 

( n-1 \ n-1 

1 ~ y ^ 9n + y ^ ^i9i 

n 

= ^707 (14) 

1 

This configuration represents the best combination 
of the original configurations that satisfies the con- 
straint of equation (3). At this point, it is ap- 
propriate to consider the significance of this con- 
straint. It constrains the scale of the final design to 
be consistent with the scales of the original configu- 
rations. Since these configurations represent small 
variations of a particular configuration type, they 
must be nearly the same size also. The constraint 
assures that the final design will also be of this scale. 
This, in turn, assures that the quantity g d - g n in 
equation (13) will be small, which is required for the 
validity of the linearity assumption. 

For many problems, the necessity of this con- 
straint is more specific. Consider, for example, the 
design of a two-dimensional supersonic inlet of the 
type shown in figure 1. The lower wall is the ba- 
sic compression surface, whereas the upper wall is 
the canopy. For this type inlet, it has been spec- 
ified that the lower surface shall consist of a dou- 
ble wedge compression ramp succeeded by a smooth 
polynomial fairing to the throat. The wedge-section 
lengths, 1 1 and I 2 , and the length of the curved fair- 
ing Z3 vary for the individual configurations but the 
overall length is constrained to be 10 units: 

Zl + Z2 + Z3 = 10 (15a) 

The canopy consists of a single wedge section of 
length Z5 with a smooth fairing of length Iq to the 
throat, with the constant 

^5 + ^6 < 6 (15b) 

Since each of the base configurations satisfies equa- 
tion (15a), the final design also satisfies this con- 
straint because of the condition specified by equa- 
tion (3). However, even though each base configura- 
tion satisfies equation (15b), equation (3) does not 
prevent the final design from violating the inequality 
constraint, since some of the coefficients may be neg- 
ative. Thus, in this method inequality constraints 
appear as soft constraints. The condition of equa- 
tion (3) only assures that the violation, if it occurs, 
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is small in magnitude. Such a result would indicate 
that a longer canopy is superior, and perhaps that 
a slight relaxing of the constraint would yield aero- 
dynamic advantages that would more than offset the 
additional weight burden. If, however, it is not fea- 
sible to lengthen the canopy, in a future calculation 
all the base configurations should have canopies with 
the optimal length of 6. 

The calculation of the coefficients a i requires per- 
haps 1 second on a modern computer, whereas the 
performance analysis calculations may require min- 
utes and perhaps hours for complex configurations or 
computer analysis codes. 

The right-hand side of equation (2) with the coef- 
ficients given by equations (4) and (12) is called the 
predicted performance. The left-hand side of equa- 
tion (2) is the desired or specified performance. The 
actual design performance is obtained by perform- 
ing the full analysis calculation for the optimum con- 
figuration, defined by equation (14). The extent to 
which this actual performance approximates p ^ de- 
pends primarily on two factors. The first of these 
is the extent to which the original set of configura- 
tions provides the types of shapes that are effective 
in approximating the required result, that is, on the 
difference between the left- and right-hand sides of 
relation (2). If the right-hand side does not closely 
approximate then additional configurations must 
be analyzed. 

The second factor regards the accuracy of the lo- 
cal linearity assumption. This condition is violated if 
the differences between the data base configurations 
are too large, if the specified design performance is 
too far beyond the performances of the analyzed con- 
figurations, or if the problem is inherently nonlinear 
even for small variations. 


Sample Problem Calculation 


An example that illustrates the previously de- 
scribed procedure is to design a wing-body combi- 
nation for a high supersonic Mach number rocket- 
powered missile. Four configurations were analyzed 
with geometries defined by the method of reference 4. 
In the equations which follow, the 11 components of 
the vectors are denoted Cj (j = 1, 2, . . . , 11) 


y(x,t) = X(x)t 

z\x , t) = A(x){[l - <r(x)]zi(t) 

+ °{x) *2 ( 0 ) 


l (- 1 

< 

t 

< 


1 V° 

< 

X 

< 

i) 


where 


zi(t) = cii\/T — t 2 (cn < 1; upper surface) 

zi(t) = —(1 — cn) V 1 — t 2 (Lower surface) 


For 


For 


For 


1*1 < c 2 < 1, 


*2 (*) = \Ai -t 2 

(Upper surface) 

*2 (*) = -\Jcl-t 2 

(Lower surface) 

C2 < |*| < 1, 


*2 (*) = 0 
a(x) = x° l 

(Both surfaces) 

|*| < c 2 , 



4 

A(x) = A(x) = 0.075 ^ hix 1 

i=l 

where the polynomial coefficients are determined 
by specifying the slope A' = ^ at three locations: 

[ c 3 (x = 0) 

X' = < C4 (x — 1) 

{ C5 (z = C 6 ) 

For 1 > |t| > C2, 

A(x) = A(x) + 0.1625 ~ C? ) 7 Z H i xi 

\ 1 “ C2 / j 

where the polynomial coefficients are determined 
as for ^ with C3, C4, C5, and eg replaced by 
C7, c 8 , eg, and cm, respectively. 

The four configurations are shown in figure 2. 
Each configuration was analyzed with the use of the 
computer code described in reference 5. The inviscid 
drag and L/D (based on inviscid drag only) were 
calculated for a = 1° and 3° at M = 3 and 4. These 
four values of L/D were taken to be the components 
of the performance vector. These vectors are shown 
in bar chart form in figure 3. 

Several properties of the results are noteworthy. 
In terms of these performance criteria, configura- 
tion 4 is significantly inferior and configuration 2 
is significantly superior to the other configurations. 
Although the difference between the performance of 
configurations 1 and 3 is small, the remaining differ- 
ences are so large that the linearity assumption is not 
expected to be a good approximation. 

Nevertheless, an improved performance vector 
was specified, and the above procedure was applied to 
design an improved configuration. Only a very slight 
improvement over the performance of configuration 2 
was specified for two reasons. First, configuration 2 is 
so clearly superior to the others that it was suspected 
to be close to optimal, and it appeared that only a 
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slight improvement could be expected. Second, the 
variation from configuration 2 was kept small to hold 
the analysis as near the linear range as possible. The 
specified values of L/D are indicated in figure 3 by 
an arrow and the letter S. 

The value indicated by the letter P in figure 3 de- 
notes the predicted performance of the improved de- 
sign, that is, the best approximation to the specified 
performance vector that can be obtained as a linear 
configuration of the four configuration performance 
vectors, subject to the constraint of equation (3). 
Since only four components are specified and four 
configuration performance vectors are available, the 
predicted performance would be exactly the same as 
the specified performance if it were not for the con- 
straint, which reduces the number of degrees of free- 
dom by one. 

The actual performance of the improved design 
indicated by the letter D is given in figure 3. If the 
calculations were perfectly linear, these values would 
coincide with the P (predicted) values. However, 
because of the large differences in the four config- 
uration performance vectors, some nonlinear effects 
exist. Consequently, the actual performance differs 
somewhat from the predicted performance and is, in 
fact, a considerable improvement over the predicted 
performance. The improved design configuration is 
shown in figure 4. As expected it varies only slightly 
from configuration 2. 


Curves of L/D versus a for the four base configu- 
rations and for the final design are shown in figure 5. 
The better configurations have curves that peak be- 
low a — 3°. Since C 2 is a nearly linear function of a 
for a > 3°, this shift in the curve maxima indicates 
a significantly lower inviscid drag (a flattening of the 
Cj) versus a relation) in the range of a = l°-2°. 

Concluding Remarks 

A method has been described for designing a con- 
figuration for a specified performance vector, based 
on results from several similar but not identical trial 
configurations, each defined by a geometry param- 
eter vector. The theory indicated that the method 
is effective provided that (1) the results for the trial 
configurations provide sufficient variation so that a 
linear combination of them approximates the speci- 
fied performance and (2) the differences between the 
performance vectors (including the specified perfor- 
mance) are sufficiently small that the linearity as- 
sumption of sensitivity analysis applies to the differ- 
ences. A computed example described the design of 
a high supersonic Mach number missile wing-body 
configuration based on results from a set of four trial 
configurations. 

NASA Langley Research Center 
Hampton, Virginia 23665-5225 
February 23, 1988 
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Appendix 

Gram-Schmidt Orthogonalization and 
Series Expansion of Difference Vectors 

Denote the norm, or length, of the vector d by |d|. 
Then, the orthonormal set is defined recursively by 
the following procedure: 


For the general term, we substitute into the summa- 
tion in equation (A3a) from equation (A4): 

i-1 j 

e i = (dj, ej) y] c jk^k 

j= 1 k=l 

i — 1 i— 1 

= d i - J2 K*> g i) 2 c jfcd* ( c jfc = 0 for k > i) 
i=i *=i 


i — 1 i-1 

= (<*«, (c jk = 0 for fc > j) 

*=lj=l 


. di 

1 1*1 

(Al) 

«2 = ^2 - {d2,ei)h 

(A2a) 

z «2 

2 N 

(A2b) 

i— 1 


= di — ) ' (di , €j ) 6j 

(A3a) 

3 = 1 


g. = iL 

1 W 

(A3b) 


i-1 

~ (AS) 

k=i 


This procedure is sufficient to derive the e* vectors 

from the d{ vectors. However, it is advantageous to where 

express the e{ vectors as linear combinations of df 


e i — c y dy 

3 = 1 


(A4) 


i-1 

Ci-I,k = g i) c ifc ( c ifc = 0 for k > j) (A6) 

3 = i 


The coefficients in equation (A4) are computed as 
follows. From equation (Al), 


Finally, using equations (AS) and (A6) with equa- 
tion (A3b) in equation (A4) yields 



From equations (A4), 



C21 


(d2,ei) 

l«2||dl| 


c ik — 



{k = i) 
(k < i ) 


(A7) 


where the coefficients are obtained recursively 
from equation (A6). Thus, with Cfc computed by 
relation (A7), equation (A4) is the result required in 
the main text (eq. (8)). 
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Figure 1. Supersonic inlet example illustrating application of constraint equation. 
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(d) Configuration 4. 
Figure 2. Concluded. 
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4 13 2 D 4 13 2 D 4132D 4132D 

Configuration Configuration Configuration Configuration 
ct= -|° a=3° a=1° a = 3° 

M = 3 M = 4 

Figure 3. Computed performance of four base configurations and design configuration, together with 
specified (S) and predicted (P) performance values. 
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